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Abstract 

The recent method of the description of classical fields in terms of Lagrange- 
Souriau form is applied to the case of source-free electromagnetic field in or- 
der to check its computational capabilities. The relevant calculations are rep- 
resented in all details and yield a useful data for comparisons of the method 
with more usual approaches in this simple and transparent case. 



In the recent papers a new general method of description of systems of classical parti- 
cles and fields has been suggested [1,2]. It seems interesting to see, however, its advantages 
(or perhaps disadvantages) in applications to concrete problems, at the first time, well 
known and easily tractable by means of more usual methods. 

As one of the simplest examples, we shah consider a source-free electromagnetic (E-M) 
field on a space-time that is admitted to be curved. The inverse infiuence of E-M field to 
the space-time is neglected and, thus, the metric is assumed to be known. Let us denote 
space-time coordinates as a:'^, electromagnetic potential as 



A. = il>adx°'. 



In the formalism under consideration field variables ipa are interpreted as quite independent 
coordinates (together with x^) of some multi-dimensional evolution space. Moreover, all 
the derivatives of field variables entering the ordinary Lagrange function describing the 
field are considered as independent coordinates as well. In our case one has the additional 
coordinates that correspond to to partial derivatives of potentials It is useful to 
introduce the following notations here: 
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In frames of formalism under consideration the physical field is described by the closed 
Lagrange- Souriau p-form a on the space with coordinates {a:'^, -^q, Xq^}. g may be locally 
represented as exact differential 

a = dO (locally), 
of the generalized Cartan form Q in turn defined by the equation 

n 
k=0 

The coefficients XAti-'-Mt-l^i ■■■«*■ have to obey certain structure equations (see [1] for 
more details) but in principle admit the following explicit local representation in terms of 
ordinary Lagrange function: 



Pfc is the permutation group of numbers 1, . . . , A; and |^^| is the signature of C is the 
ordinary Lagrange function. 

For the vacuum electromagnetic field it is convenient to choose the Lagrange function 
in the form 



(if Xafj. are interpreted as it will be reduced to the usual representation) and we have 
respectively 



4 1 

=^i^ri^2i^si^4 J2 ^^ai A • • ■ A 6iP^, A dx'^'^^ A---Adx 

k=0 

=<^MiM2M3M4(l ■ Ldx'^' A dx'^^ A dx'^'' A dx'^'' 
4 I 

— • L'^il"! 6^^^ A dxf"^ A dxf"^ A dx''* 
where we have denoted for convenience 



fJ-4 



Vi^ii^2 ^1^1 1-^2 i^ai^A^"^^ A dx^ , 
Vfil — ^^1/^2/^3/^4 ^'^'^^ ^ dx^^ A dx^'^ — dx^^ A ^^^^2' 
77 = e^^^^^^^^dxf^^ A dxf^^ A dxf^^ A dx'^'' = dxf^^ A r?^^. 
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Other Us vanish. One has further 

L = C, 



~ ~ ^XaMX/3'/V 99 9 



7-MiM2|aia2 _ 1 
^ ~ 2 
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1 / ^M2|a2 ^Ml|a2 

aiM2 



Meanwhile 



^M2|«2 = ^ 4^^^^^«2[/3^.]M2 =4y3^^«2[«l^Ml]M2 



and therefore 



_^^«iMi^«2M2 _ ^g'^^l'^gl'^''^ 



Thus 



It is essential to mention that this generalized Cartan form contains, besides antisym- 
metric part of Xo/x, its symmetric part as well. 

The Euler-Lagrange equations in frames of the formahsm under consideration take 
the form 

^^*{X^Ja) = 

which equation has to be observed for arbitrary vector field X on the evolution space. Here 
denotes its lift to the jet-1 bundle, is the pull-back of the immersion $^ of the 
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space-time into the jet-1 bundle which is described in the coordinates {x^^, ijja, Xa/j.} by the 
equations 



_ 

for some functions ^a{x). The latter determines the solutions of the Euler-Lagrange 
equations if and only if after such a 'substitution' the TV-form (TV is the dimension of 
space-time) J a vanishes identically for every X (the pull-back of the differential form 
is locally nothing else but the result of substitutions of the representations of coordinates 
of the space 'to where' in terms of coordinates 'from where'). 

Now we shall calculate the lift of the vector field X to the jet-1 bundle following the 
work [6]. It can be done as follows. 

The vector X can be represented in the form 

^ = n^)i + Ca(^,^)-^ 



dx>^ dlpa 

(^ and ( are arbitrary functions). This vector field is not of the most general form since 
^'s do not depend on ipa- But this is necessary restriction because otherwise infinitesimal 
transformation of evolution space produced by X (see below) would introduce into the 
functions on space-time {gafs in our case) explicit dependence on field variables il>a. 

The vector field X on the evolution space is associated with the infinitesimal motions 
of the space in accordance with equations 



X 



+ eC'ix) + 0(e) = £^ + 0(e), 



i>a ^i>c + eCc{x,i') + o{e) ='0^ + o(e). 
If XaiJ. were interpreted by means of equation 

then the 'lifted' transformed Xa^ {i-^- Xa^) were defined by means of a similar equation 

#a = Xaf^dx'^. 

Using the explicit form of the expressions for x^^, tp^ one can find the following trans- 
formation law of Xan'- 

The complete transformation expressed above in terms of jet bundle coordinates 
{x^,ipa,XaiJ.} corresponds to the vector 
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which is just the hft of the vector X to the jet-1 bundle. It turns out to be convenient 
to represent it as a sum of two 'independent' vector fields 



It is clear that Euler-Lagrange equations must be obeyed with ether X^ = Xj or X^ = X^ 
choices independently. 

Now we need calculate a = dO using the explicit representation of given above. The 
following obvious equations are useful: 



dSipa = -dxcp. A dx^^, dr) = dr]^ = dr]^^^^ = 0. 



de= 4 rfXap A 77 

+ 16 v^5"iI/35-lMi dxp. A SiP^, A rjf., 

- 16 v^5"^'^5"'''^ Xp. rfX«,Mo A dx''^ A 77^, 

+ 6{v^5"^''^5''^"^},Mo A SiP^, A dx''' A Vf., 



Some simplifications can be achieved by means of equations 



X' = XI + X^ 



c 



where 




One has 



dxf"' A 77^, 

dX^ A 77yL(,iyU2 



They yield 



-4{v^5"^'V'^},;.X/3.^^«. AT? 
+ 16 v^5"^'^5"'''^ A A 77^, 
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+ 8 V^^-iI/^^/^^]"^ d;^^^^ A A 

+ 16 v^5"'^5"'^ A .5^^ A Vi. 
+ 4{v^/['^5^1'^},. 5^/3 A Si>^ A 77;, 
+ 8 v^/'"5'^l" A A 7?^. 

This is an explicit representation of Lagrange-Sonrian form for electromagnetic field. 
As well as generalized Cartan form it involves symmetric part of Xa/j. because entering the 
last term dx(pu) certainly cannot be cancelled out. 

Now we can calculate the contractions ZjJ a and Z^A a. It must be emphasized how- 
ever that these contractions produce the Euler-Lagrange equations when they are restricted 
to the image of the space-time under the immersion The latter however possesses the 
property to annihilate Sipa- (Sipa) = 0. Thus, in order to construct the Euler-Lagrange 
equations it is sufficient to keep only those terms of the contractions that do not contain 
the factor Sipa- But all the terms of d9 representation do contain it at least once. Hence, 
the vectors or Z^ have to be contracted just with such a factor and the term may not 
contain other factors Sipa if it contributes the resulting equations. 

One has obviously 

Zl_\dx»/^ 

z^Jsi,^ 

Similarly 

Z^J Sipa = Ccix,^)) 
Z\Ay]^ = ^. 

It follows also from these equations and the above remark that in our case 

^^\z\Ag) = ^ =4> ^^\z\Ag) = ^. 
and therefore it is sufficient to calculate Z\Ag mod Sipa- 
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The contraction with Lagrange-Souriau form a = dO will be the following: 

- Ca 16 V^g'^^^g"^'' dx(3. A 
It remains to substitute xpty = di,^p{x) to obtain the following version 

+4v^5"'V''^rf^.*/3A77;, 
+2 v^/[-5^1" (^oi,*^ A r;^ = 

of the vacuum Maxwell equation. 

Is it a true one? It is not so obvious at the moment. 

Let us note however that the last term of the above equation vanishes identically 
because 

dd„^f3 AVi^ = dpu'^p dxP AVf. = ldf,^<fp r] 

is symmetric in fj,,u and the factor gf^l'^g^^]" is antisymmetric. Then the implied equation 
can be expressed in term of antisymmetric E-M stress tensor 

alone simultaneously removing antisymmetrizations from the metric tensor aggregates: 

{ V^g'^^g"^ ],^Fp,rj + A V^gg'^Pg"^ dFp, A 77^ = ^ 
{ V^5"V^},M Fp, + V^gg^'^'g'"' d^Fp, = 
ga, { V^gg^^g""^ } f^^ + ^gg^^^ d^F^, = 

gc.AV^9''^9'"'Fp.],^ 

where we have defined as usually F°''^ = g^'^g'^'^Fp^. This is just the vacuum Maxwell 
equations 

Vf.F'^'^ = 

that follows from the wellknown definition of covariant derivative: 

a^F"^ = Vf.F'"' - VJ^pFP"" - V'^pFP'^ =^ 
^M^"^ = V^F"^ - T'/^pFP'^, 
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We have obtained Maxwell equations in 'apparently covariant' (with respect to co- 
ordinate transformations) form. The general definition of Q can be seen as 'apparently 
covariant' as well. Meanwhile the above expression of the Lagrange-Souriau form obtained 
by straightforward calculation does not reveal such a feature. It makes sense to obtain an 
'apparently covariant' representation of a also. 

Let us transform a 



eliminating derivatives of metric tensor by means of formulae 



One has 



similarly 



Pill 
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a/3 _ _-na „p/3 _ "n/B pa 
1 HP J HP " 



MP- 



:^(^_^«P^M[-r^]^_^PM^«[/3r- 



and therefore 



8v^ 
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Next we introduce "formally covariant" objects 

and their "formally covariant differential" 

D X(Ti^ dX(Ti^ Xpij^ dx^ X^y^u Tyr^p dx^ . 
One has therefore 

dXoM = dX^^ + di^p + i^p dT^^, 
where it is necessary to substitute for 

di>p = 8tpp + Xf3p dxP = 8tpp + Xpp dxP + -^^r^p dx^ . 
It is worthwhile to mention that 

X\ap] = ^[a/x], dxia/j,] = dX^a^] = + Xia^]T'^^pdxP + X[p^]T'^pdxP 

and one has for 

=2g'^f'g'^^DX[p,] A /\rjp- i5"^5"^^[/3.]r;:/^« A TJ 
=2g-^g'^^DX^p,^ A A - i^^^^^^Xf^^] r^^5^« A 7? 

=2 5"^5^'^Z)X[^,] A A 7?^ - 5"'^5"''^^[/3.]r^^5^« A V. 
Meanwhile without antisymmetrization one has 

dxp. = DXp, + T}^ 8i^^ + {^^ d,Tl^ + x^p T}^ + r^pr;^, + x,, v}^ + r^p} dx'^ 

and therefore 

A A ={i:) + 5^^} A A r]f, 

- {^7 ^pr^^ + x,p rj^ + r^/j, + x,. r^^ + Xp„ vz^} x 

^-^a A dxP A 7?;, 

X^, + 5^^} A A 77^ 

- i{^7 ^/.r^, + r]^r^, + 2 x^^^ r^^^ + x^. r-^}5^« a 77, 

rf^^^ A 8^^ A 77;, =/[^5'^'" Xp^ A a 77;, + /'^^'^'^r^^ 8^^ A A 77^, 
- i/^5'^"{^7 ^[,r3]^ + TlTl^^}8i^^ A 7?. 
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The latter expression contains precisely the Riemann tensor 

and, after contraction with metric, the Ricci tensor = W^^°'^ enter it as follows 

/I-^/^]" A A 77^ =/[-5'^l" X^, A A 77^ + /["^'''"r^^ A A 77^ 

+ lij^WStjj^ A 7?. 

As a result one has the following expression of the the Lagrange- Souriau form 

- g^i'-g^'^rg/^^ A A 7?^ 

+ 2 ^[/3.] A A 77^ - g"'V^r-,X^,^^^ A77 

+ /l^^'^^" D Xp^ A A 7?^ + g^^^g^^^^TlJ^^^8iJ^^r)^ 
+ li>^R^''Si>^ A 7?. 

Underlined terms are mutually cancelled and 

i^^i?'>"^^^« A77 
+ 2 ^"^^^'^ A A 7?^ + /I-^'^]" A A 7?^ 

Let us consider another method of solution of the same task and begin with the above 
representation of the generalized Cartan form 6. Due to the obvious indices symmetry 
reasons it can be recast in 'apparently covariant' form as follows: 

+ eV^g'^'^'g^''^' Sip^, A Sip^, A 7?^,^,. 

The exterior differential d of this expression whose all the indices are contracted can 
be replaced by the 'covariant' differential D which is constructed in accordance with index 
structures of object in a way usual for tensor calculus and commutes with multiplication 
to the metric tensor: 
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dB = De = 

One has obviously 



HP 



AT] = 



|(\/^),[mi'7m2] - i 



^ P[M1 'M2] 



Further, 
and therefore 
But 
and 



DSip^ = DDip„ - DX^f, A dx^ - X^^^Ddx^. 



Ddx'' = ddxf" + V^JxP Adx" = 



DS^a = DD^a - DX^f, A dx^. 
We have the following representation: 

d0 = 2v^5«^5^'^ + i^X,^^]) A V 

+ 16v^5"i^5^'^i {dX^p^] a 6iP^, + DSiP^,') A r]^, 
+ 127^5"^'^^ 5^^^"^ D8tP^, A SiP^, A 77^,^,. 

It contains the terms which are not 'apparently covariant' and have to be expanded: 

DSi>ai A ?7^i = DDtp^^ A T]f,, - |M«j^j a 77, 

DSipa^ A Sipa^ A 

= DD^a^ A A 77^j^2 + ^^«M3 A A dx'''' A r?^^^^ 

= DDipa^ A Sipa^ A r]p,^p,^ + iDX«j^j A Sip^^ A r?^^ - ^DX^^^,^ A Sip^^ A r]^,^. 
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This yields 

de = 4v^5"^5''''X[^,.] DXi^ij,] A r] from the first line 

- A^/^g''^^ g"'^'X[i3,j]DX,yf^ A r) from the second fine 

+ 16\/^^g"'^g'^^DX^pi,] A Sipa A 77^ from the second line 
+ 'iV^g"'''' g^^'^^ DX^^f,^ A Sij)^^ A Tj^^ from the third line 

- 4v^5"'^'5^'"' ^^aiM2 A A 7?^, from the third line 
+ \Q^/^g"^ g"^X\^i3j^]^DD'^a A rj^ from the second line 

+ 12v^5"i^i 5'''"' i^i^^ai A /\ ^MiM2 from the third line 
= le^gg^^g^^'^DXy,,,^ A A 7?^ 

+ 8v^5^'"5^'" A A 7?^ 

+ lev^^-^^^^'^Xi^,]^)^)^^ A 7?^ 

+ 12 v^^"!'^!^'^^"^ Z>£»^«, A A 7?^,^,. 

The duplicated covariant differentials can be expressed in terms of Riemann tensor as 
follows 

=d{#« - i>pTi^dx^} - ri^dx^ A {di,p - i>^r}^dxp} 

= - i'pdr'l^ A dxf" - Ti^dtPp A dxf" - Ti^dx'' A + ^^F^^r^^da:'^ A dx^ 

= - i'p{dri, + n^r^.^dx^} A dx^ = -i>f3{d^r^^^^ + r^i/ij^}^^" ^ 
= - I'^nRip^dx" A dx^ 

and since 

dx^ A dx^ A 7?^ = 0, dx^ A dx'^ A 77^^^, = l^f^/J^^j^? 

= 16V^g''"g"''DX[p„] A .5^^ A 7?^ 
+ 8v^/[^5'^l" A A 7?^ 

= 16v^5"''5"'^i^^[/3.] A 6iP^ A 7?^ 
+ 8v^/[^5'^l" A A 7?^ 

+ v^i?^"^^5^« A 77 

that coincides with previously obtained representation. 

What comments would be relevant to the calculations exposed above? One may 
suggest a number of them and all of a negative sort. 
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At first, the replacing of tlie potential 1-form A. = ij^adx"" being the invariant geo- 
metric object in favor of the indexed components -^q,, Xa^ seems to be disadvantage badly 
conformed with the current tendency to work in invariant notations. Besides, this makes 
much more complicated the index structure of the expressions under consideration. 

The second disadvantage related to the former one is a too involved index structure 
of equations comparatively with the problem considered. It is of course more technical 
than principal feature but in practice it may became a serious obstruction on the way of 
application of method to real and much more involved than this example problems. From 
the point of view of computational convenience the methods seems to be not attractive. 
In particular, it seems doubtful to recommend implementation of the method to the real 
computer algebra systems. In fact, if such a system is assumed to be able to work out 
with arbitrary Lagrangians it should be equipped with "absolute" capabilities to process 
numerous contractions of indices, taking into account and properly exploiting any type of 
relevant symmetries etc. 

Next, it seems rather unattractive that the basic object o includes effectively the terms 
that are related to symmetric part of the partial derivatives of E-M potential. Perhaps, it 
is connected with some more deep physical background, but there is no clear evidence of 
this yet. 

The last relevant remark is the following: one can see that a large part of calculations 
leading to the Lagrange-Souriau form a is in fact useless since after 'contractions' and 
'restrictions to space-time manifold' a number of terms required some efforts to calculate 
them are merely dropped. 

One has to emphasize here however that the methods of descriptions of physical sys- 
tems in term of evolution spaces and generalizations of Cartan form are destined mostly 
to study their symmetries rather than to derive evolution equations [3-6]. Thus it is rea- 
sonable to outline how the symmetries of E-M field considered as dynamical system are 
described in terms of Lagrange-Souriau form formalism. 

The 'apparently covariant' representation of a obtained above immediately implies 
its invariance with respect to coordinate transformations '^(rc), provided i\\dX 

i'a, Xa/j. are transformed by the tensorial rules in accordance with standard interpretations 
of their indices. It is not necessary however to restrict the class of admissible transforma- 
tions of il>a, XaiJ. in such a way and much more general transformations may be involved in 
the play. Then more general group of symmetries of the system might be found. We shall 
consider here less general problem and restrict ourselves to the case when the coordinate 
part of transformation is trivial, i.e. x^{x) = x^. We are guaranteed that such symmetry 
subgroup is still non-trivial since it contains the gauge transformation of E-M potential at 
least. Thus we shall consider the transformation of the Lagrange-Souriau form 

+ 8 v^{2 5"[^5-l^ + /V'" }dXp. A Si,^ A 7?^ 

under the map 
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X' 



i>c. ^ i>c. = i>c. + i^c.{x, i>, x), x«M ^ x«m = x«m + Xc^ai^, i>, x)- 



for arbitrary functions -0, x Tlie ordinary E-M gauge transformation corresponds to ip„ = 

^gjf , XaiJ. = dx^ik'l ^'^^ arbitrary function A{x). It is easy to see that a is invariant with 
respect to it implying that gauge transformation is a symmetry. One has obviously 



- a = - 4{ v^5"[^5-l^},^ Xf3. <5^« + Xf3. + Xp. St 



A 7] 



where the evidently generalized notation is used: dip^ = dip^ — Xafj.^^^^- The components 
ip, X have to be interpreted here as functions of x, ip, x- 

Now it must be noted that Noethern symmetries are defined in [1,2] just as dilfeomor- 
phisms of evolution space leaving a unaltered, i.e. as solutions ^ of the equation 



^ a = a. 



Strictly speaking, this restriction is too strong and it would be sufficient for the form 
(7 — (j to vanish after contraction with the hft of arbitrary vector field being restricted 
to the image of space-time under the immersion $ and on solutions of the Euler-Lagrange 
equations only in order to map a solution to another solution that is required from a 
symmetry map. On the other hand, if a is invariant itself this implies, at list locally, 
the relation ^ Q — Q = dVl ior some and then the Noethern current can be easily 
constructed. Thus even if we adopt the above more restrictive definition of the symmetry, 
the corresponding equations in explicit form will be rather complicated nonlinear system 
of partial derivatives equations. Its only obvious general property is complete integrability 
(both a and a are automatically closed) that means here that no additional information 
can be obtained with the help of exterior differentiation of the equation. Of course, when 
the problem is considered in the whole generality and arbitrary diffeomorphisms of the 
space-time x^ ^ x^ = x^{x), are admitted the equations become much more involved. 
Certainly, one may obtain some their partial solutions but to describe the whole set of 
symmetries seems to be unreahstic goal, at least if the above form of the equations is used. 

In conclusion it is worthwhile to notice that the point of view followed in these notes 
concerns mostly with the computational efficiency of the Lagrange- Souriau form method. 
It is quite possible that its other kind features would reveal some important advantages of 
the method as well. 
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Appendix 

The shortest way from Lagrangian to Maxwell equations seems to be the following: 
we may define L[A\ = ■^(dA) A *dA. where * is Hodge dualizing and then 

6L[A, SA] = ^^1-^ + ^^-^J = lidSA) A *dA+ lid A) A *dSA 

de |e=o 

= ^SAa d*dA+ ld{SAA *dA) + li*dA) A dSA 
= SAAd* dA + d{SAA *dA). 

The factor d * dA coupled with SA in the first term is the left side of Maxwell equations, 
SAA *dA is the boundary term. 
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